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Abstract 

We show that the model which naturally explains the u 7^ d asymmetry in the 
nucleon and is in quantitative agreement with the Gottfried sum rule data, 
also predicts that in the proton Am > > As > Ad and Au — Ad > d — u > 0. 
At the input scale, these results can be derived even analytically. Thus the 
violation of the flavor symmetry is more serious in the polarized case than in 
the unpolarized case. In contrast, many recent analyses of the polarized data 
have made a simplifying assumption that all the three Ag's have the same 
sign and magnitude. We point out the need to redo these analyses, allowing 
for the alternate scenario as described above. We present predictions of the 
model for the W~ asymmetry in polarized pp scattering, which can be tested 
at RHIC; these are quite different from those available in the literature. 
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Several comprehensive analyses of the polarized deep inelastic scattering (DIS) data, 
based on next-to- leading-order quantum chromodynamics (QCD), have appeared recently 
P|~P|. In these analyses the polarized parton density functions (PDFs) are either written 



in terms of the well-known parameterizations of the unpolarized PDFs or parameterized 
independently, and the unknown parameters are determined by fitting the polarized DIS 
data. Additional simplifying assumptions are often made; the one that has been widely 
used in the literature fl|-|§ is 

Au = Ad = XAs, (1) 

with a positive A which is usually set equal to unity. Recently the HERMES and SMC 
collaborations |T2]JT3[| too analyzed their inclusive and semi-inclusive DIS data assuming all 



Ag's to be of the same sign. The same assumption has also been used to make predictions 
for future accelerators; see e.g. [Q. In some analyses [pj] , the nonsinglet PDFs Aq% and 



Aq 8 are assumed to differ only by a constant multiplicative factor (see (21) below). 

In this paper, we examine these simplifying assumptions made in the literature. This 
is important because a similar ad hoc assumption about the flavor decomposition of the 
unpolarized antiquark sea, u = d, turned out to be wrong when accurate data on muon DIS 
became available ||15|| , and the global analyses of the unpolarized data had to be redone. 



Here we derive a series of inequalities satisfied by the PDFs and point out the need to redo 
the global analyses of the polarized DIS data in the light of these inequalities, allowing 
in particular for the violation of the flavor symmetry in the polarized antiquark sea; see 
(20) below. We then present predictions of our model, which can be tested in polarized pp 
scattering at RHIC, BNL. Finally, we describe other recent works on flavor asymmetry of 
polarized sea distributions. 

We use the framework of the statistical model for polarized and unpolarized struc- 
ture functions and PDFs of the proton and the neutron, which was presented recently 
Ttj|jrT| ]. This model provided a natural explanation of the u ^ d asymmetry in the nu- 



cleon and was in quantitative agreement with the Gottfried sum rule data. Additionally, 



it reproduced the data on F 2 p (x, Q 2 ) for 0.00001 < x < 1 and 2.5 < Q 2 < 5000 GeV 2 , 
F$(x) — F£(x), F^{x) / F^x), xg(x), d(x) — u(x), d(x)/u(x), the fractional momentum of 
charged partons and the polarized structure functions g{' n (x), at various Q 2 . Out of those, 
only the F$ and (F$ — F£) data, both at Q 2 = 4 GeV 2 only, were used as an input to fix 
the model parameters, and all other results served as model predictions. In particular, the 
d(x)/u(x) ratio in the limit x — > 1 turned out to be 0.22 in good agreement with the QCD 
prediction 0.2 |18[| . At the input scale (Q 2 = Ql = M 2 , where M is the nucleon mass), 



all xq(x) and xq(x) distributions were found to be valence- like, and xg(x) was found to be 
constant in the limit x — > 0. Thus the total number of gluons was logarithmically divergent 
providing a strong a posteriori justification for the statistical model ansatz ||17|| . Contrary 
to the common practice, the polarized and unpolarized data were reproduced in a single 
framework and the simplifying assumption of charge symmetry was not made. Here we 
further explore the predictive power of the model. 

If 7"Wa)T(J.) denotes the number of quarks (antiquarks) of flavor a and spin parallel (an- 
tiparallel) to the proton spin, then any model of PDFs in the proton has to satisfy the 
following constraints: 

(2) 
(3) 
(4) 
(5) 
(6) 
(7) 

The RHSs of (5)-(7) have been measured by several groups. We use (Am + Au) = 0.83 ± 
0.03, (Ad + Ad) = -0.43 ± 0.03, (As + As) = -0.10 ± 0.03; see |TJ|. The parton numbers 
n a(a)i(i) in (2)-(7) are obtained by integrating the appropriate number density dn/dx over 
x. The various A's are also ^-integrated quantities. 
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The RHSs of (2)-(4) are clearly Q 2 -independent. The RHSs of (5)-(7) are also Q 2 - 
independent in the jet and Adler-Bardeen (AB) schemes: Recall that the nonsinglets Ag 3 = 
(Au + Au)-(Ad + Ad) and Ag 8 = {Au + Au) + [Ad + Ad) -2(As + As) are (^-independent 
in all renormalization schemes because of the conservation of the nonsinglet axial vector 
current, and the singlet AS = (Am + Au) + (Ad + Ad) + (As + As) is Q 2 -independent in 
the jet and AB schemes because of the Adler-Bardeen theorem |[20| . As a result, (Am + 
Am), (Ad + Ad) and (As + As) which can be expressed as linear combinations of Ag3, Ag§ 
and AS, are also (^-independent in these two schemes. In the MS scheme, on the other 
hand, AS is Q 2 -independent at the leading order and only weakly Q 2 -dependent at the 
next-to-leading order. Empirically too AS is found to be almost Q ^independent; see e.g. 



Fig. 5 of PI| . Hence in the MS scheme the RHSs of (5)- (7) are expected to be nearly 
Q 2 -independent. 

We now show how the statistical model naturally leads to a violation of the flavor symme- 
try in the unpolarized and polarized seas in the nucleon. Consider the following 6 equations: 



2n uT - 2n ul = 2.83, 


(8) 


2n ul - 2n m = 1.17, 


(9) 


2n d] - 2n Sl = 0.57, 


(10) 


2n dl - 2n d - T = 1.43, 


(11) 


2n s1 - 2n n = -0.10, 


(12) 


2n sl -2n 5T = 0.10. 


(13) 



These are obtained from (2)-(7) by linearly combining the latter set of equations in pairs. 
E.g. (8) and (9) are obtained by adding or subtracting (2) and (5). 



It was shown in |T6| that the parton number density dn/dx in the infinite-momentum 
frame, at the input scale, is given by 

dn _ M 2 x f M / 2 dE dn 

where 



dn/dE = g f(E) (VE 2 /2n 2 + aR 2 E + bR), (14b) 

is the density in the nucleon rest frame. Here M is the nucleon mass, E is the parton energy 
in the nucleon rest frame, g is the spin-color degeneracy factor, f(E) is the usual Fermi or 
Bose distribution function f(E) = {exp[(E — fi) /T] ±1} -1 , V is the nucleon volume and R is 
the radius of a sphere with volume V. The three terms in (14b) are the volume, surface and 
curvature terms, respectively; in the thermodynamic limit only the first survives. The two 
free parameters a and b in (14b) were determined in [|17] by fitting the structure function 



F2(x,Q ) data at Q = 4 GeV . Their values as well as the values of the temperature (T) 



and chemical potential (/i) which get determined due to (2)- (7), were given in |T7 
At the input scale, with the help of (14), (8) can be written in a full form as 



r i M 2 r r M / 2 dE 

/ dx ^ ^ g (VE 2 /2tt 2 + aR 2 E + bR) 

JO 2 JxM/2 E 2 



2 2 



a j3{E-n u] ) _|_ j e P(E-fj, n ) _|_ l 



2.83 . 
(15) 



It is straightforward to show that the chemical potentials for quarks and antiquarks satisfy 
the relations 

/^t = ( 16a ) 
Ni = ( 16b ) 

So it follows from (15) and (16b) that fi u ^ > 0. Similar arguments show that /i^, /i^ 
and n s i are positive and /z^ is negative. Moreover, since the RHSs of (12) and (13) differ 
only in sign, we have = —fi s i- Since RHSs of (8)-(13) can be arranged as 2.83 > 1.43 > 
1.17 > 0.57 > 0.10 > —0.10, the corresponding chemical potentials satisfy 

/M > m > /M > M > (AM = AM) > > (AM = AM) > AM > AM > > AM- ( 17 ) 



It will be useful to recall the actual values of the /i's given in |L7j- They are (in MeV) 
fjL u 1 = 210, = 106, fi u i = 86, /vr = 42, = 7, /x s f = —7. /i's for the antiquarks 
follow from (16). [The RHSs of (8)-(13) are sufficiently different from each other so that the 



experimental errors in (Aq + Aq), quoted above, will not alter the ordering in (17).] (17) 
together with (14) yields, at the input scale Q\ (= M 2 = 0.88 GeV 2 ): 

n«t > n d[ > n ul > n d] > (n sJ . = n H ) > (n sT = n^) > n Sl > > > n U [ > 0. (18) 

As a check, it is easy to verify that (18) reproduces the correct signs of the RHSs of (2)-(7). 
Notice the symmetric arrangement of the /i's in (17) and the consequent arrangement of the 
n's in (18). 

To recapitulate, the statistical model provides a quantitative method to incorporate the 
effects of the Pauli exclusion principle into the PDFs: the RHSs of the number constraints 
(2)-(7) or equivalently (8)-(13), force the various chemical potentials and hence the parton 
distributions to be arranged as in (17) and (18), respectively, at the input scale. 

Further consequences of (18) are easy to derive: (Note n q = n q ^+n q i and Aq = n q ^—n q i.) 

(a) The general positivity constraints on the polarized and unpolarized PDFs: \Aq\ < n q 
are satisfied trivially. 

(b) Am > 0, Ad < 0, As < 0. 

(c) Am > 0, Ad < 0, As < 0. This is in contrast to the assumption (1) made in the 
literature [|l]^|r^-[rj| that all the three Ag's have the same sign. 



(d) Au v = Au — Au = n u ^ — n u i — + n^i > 0, because the two n u terms are too small 
compared to the two n u terms (see (18)) to change the sign of the RHS. 

(e) Ad v = Ad — Ad = — — + < 0, because the two nj terms are too small 
compared to the two terms (see (18)) to change the sign of the RHS. 

(f) As,, = As - As = 0. 

(g) Ag 3 = (Au + Au) - (Ad + Ad) > 0; see (b)-(c). 

(h) ng > n a which leads to the Gottfried sum rule violation. Thus the statistical model 



naturally leads to the u ^ d asymmetry in the unpolarized sea ]IB| . Moreover, it was shown 
in ]T7|] that the model is in quantitative agreement with the data on (Ff — Fg) vs x and the 
Gottfried sum Sq. 
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(i) Am — Ad > ng — n„ > 0. Thus the violation of the flavor symmetry is more serious in 
the polarized case than in the unpolarized case. 

(j) Ad — As = na\ — — n s -\ + n s ^ < 0, because n s <[ and n s \ tend to cancel each other, 
unlike and n^. Combining this result with (b) above, one gets \ Ad\ > \ As\, and 

Aw > > As > Ad. (19) 

(k) Ad — As = njj — nj| — + n^i < 0, because and n S | tend to cancel each other, 
unlike and n^. Combining this result with (c) above, one gets \Ad\ > \ As\, and 

Am > > As > Ad. (20) 

We have derived the results (a)-(k) analytically, at the input scale. They are borne out by 
actual numerical calculations; see Fig. 1 which shows our polarized PDFs at the input scale 
Qq = M 2 = 0.88 GeV 2 . We have evolved our polarized PDFs in the next-to-leading-order 



QCD, in the MS scheme, in the range Ql < Q 2 < 6500 GeV 2 . We find that the results (a)- 
(k) are valid throughout this range. Figure 2 shows that the violation of the flavor symmetry 
is more serious in the polarized case than in the unpolarized case, throughout this range. 
Incidentally, we have examined another simplifying assumption made e.g. in ]TT], namely 



Aq 3 (x,Q 2 ) = C Aq 8 (x,Q 2 ), (21) 

where C is a constant independent of x and Q 2 . The present model predicts that (21) is not 
justified (Fig. 1). 

The statistical model makes concrete predictions for various asymmetries in polarized pp 
scattering, which can be tested at RHIC. For example, parity-violating single- and double- 
spin asymmetries for W production in the reactions jfp — > W ± X and — > W^X 
respectively, are given by p2| . [23[| 

PV + = Au{x a , Myy) d(x b , Myy) - Ad(x a) M ^ ) u(x b , M 2 ,) 

L 1 ' uix^M^dix^M^ + dix^M^uix^M 2 ,) ' 1 ' 

A P L V (W-) = - A * d + Ad * (23) 
L v ; ud+du ' v ; 



.PV/rr - Am d — u Ad — Ad u + d Au 



A PV (W + ) = ^ i ^ , s 

LLK ' u d- Au Ad + du- Ad Au y ' 

a py , . u Ad — Au d — d Au + Ad u 

All [W } = ud-AuAd + du-Ad AS' ( } 

where x a = y/re y , Xb = \fre~ y , r = M^/s, y is the rapidity of W and >/s is the pp center- 
of-mass energy. The arguments x a , Xb and Mjy are suppressed in (23)-(25) for brevity of 
notation. 

In the present model, Am and Am are positive and Ad and Ad are negative (see (b), (c) 
and Fig. 1). Also note that Am < u and \Ad\ < d (see (a) above). Hence it is straight- 
forward to show that < A^ V (W + ) < 1. Similarly, -1 < A^ V (W-) < 0, Af%(W + ) > 
0, AlY{W~) < 0. It is somewhat tedious but again straightforward to show using (22)-(25) 

that A£¥(W+) > A P L V {W + ) and \A P ¥( W ~)\ > \ a l V ( w ~)\- A q uick and crude wa Y to 
convince oneself that A P j j (W + ) > A PV (W + ) is to ignore the (small) "AA" terms in the 
denominator of (24), which makes the denominators of (22) and (24) identical, and then to 
compare their numerators. In fact, at y — (or x a = Xb), A P Y(W + ) is seen to be almost 
twice as big as A^ V (W + ). 

Figure 3 shows our predictions for A pv and A^Y for W~ production in polarized pp 
scattering at ^/s = 500 GeV, as a function of the rapidity y. The above inequalities for 
A PV (W~) and A P Y(W~), which we derived analytically here are borne out by the actual 
numerical results in Fig. 3. Also shown for comparison are results reported in |[23|| . These 



are based on the parameterizations of polarized PDFs given in PJ3"U24|. Asymmetries for W 
production are sensitive to the sign of Am which is positive in the present model, negative 



in P,|24| and x-dependent in 0. The recent work of de Florian and Sassot |25j has yielded 
a clear preference for a positive Am distribution. 



As stated earlier, the HERMES and SMC collaborations ||T2"| , |i"3"f analyzed their inclusive 
and semi-inclusive DIS data assuming all Ag's to be of the same sign. Recently, Morii and 
Yamanishi have reanalyzed these data and have estimated Ad(x) — Au(x) at Q 2 = 4 



GeV 2 . It is evident from their Fig. 1 that Au(x) — Ad(x) is positive and has a peak at 
x ~ 0.06 where x(Au(x) — Ad(x)) is ~ 0.05. All these observations are consistent with our 
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Fig. 2. 

Another model which is able to generate flavor asymmetric polarized antiquark sea is the 
chiral quark soliton model (CQSM) p7|-p2|. Results in our Fig. 2 are strikingly similar to 



those in . This is remarkable because the physics inputs of the two models are quite 

different. It is also noteworthy that the origin of the u 7^ d and Am 7^ Ad asymmetries is quite 
simple in the statistical model. While the role of gluons is yet to be understood in CQSM, 
the statistical model predicts a positive Ag(x, Q 2 ). The pion cloud model also gives rise to 
the u 7^ d asymmetry, [for a recent review, see [p3| 1, and there have been some attempts to 
generate polarization by including spin-1 resonances in that model. These attempts have 
been commented upon in |3yj31| . Recently, Gliick and Reya p4| have discussed the issue of 



flavor asymmetry, in a phenomenological way making use of the Pauli exclusion principle. 
We recall that the statistical model [|T6| , p~T|| provides a quantitative method to incorporate 
the effects of the Pauli exclusion principle into the PDFs. 

We have treated all partons as massless: m u = = m s = 0. If m s is taken to be 
nonzero, then (14) will have to be generalized, but the parton densities still have to satisfy 
(2)- (7) and equivalently (8)- (13). So it is not obvious how that will affect the symmetric 
arrangement of the /x's in (17) and the consequent arrangement of the n's in (18), at the 
input scale. This is a nontrivial problem which needs to be investigated further. 

In conclusion, we have derived, on rather general grounds, a series of inequalities for the 
polarized PDFs; see (a)-(k) above. This points to the need to redo the analyses [p]^8], p!2|JT3|| 
of polarized data, allowing for the alternate scenario as in (19)-(20). Some of the inequalities 
can be tested in the forthcoming spin-physics program at RHIC, BNL. To illustrate, we have 
given our predictions for the W~ asymmetries; these are quite different from those available 
in the literature. 

I would like to acknowledge the hospitality of the Nuclear Theory Center, Indiana Uni- 
versity where this work was initiated. 
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FIG. 1. Polarized PDFs at Q 2 = Q 2 (= M 2 = 0.88 GeV 2 ). 
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FIG. 3. Predictions for asymmetries in W~ production in polarized pp scattering at y/s = 500 
GeV vs rapidity y. Solid curves: present model. Other curves: Long-dashed [24], short-dashed [2], 
dotted and dot-dashed [3] for two separate parameterizations — as reported in [23]. 
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